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I A connection is establisehd between the results of a potential quark model with dynamical masses 

;h ! 

Ph' and the l/Nc expansion mass formula used in the description of baryon resonances. It is shown 



that a remarkable compatibility exists between the two methods. 
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I. INTRODUCTION 



So far the standard approach to baryon spectroscopy is the constituent quark model 
where the Hamiltonian contains a spin independent part formed of the kinetic plus the 
confinement energies and a spin dependent part given by a hyperfine interaction. The latter 
can be either due to one gluon exchange or to Goldstone boson exchange between quarks, 
or it can be an instanton induced interaction. The results are naturally model dependent. 

It is therefore very important to develop model independent methods that can help in 
alternatively understanding baryon spectroscopy and that can support quark model assump- 
tions. Large Nc QCD offers such a method. In 1974 't Hooft proposed to generalize QCD 
from SU(3) to SU(A^c) where Nc is an arbitrary number of colors and suggested a per- 
turbative expansion in the parameter 1/Nc, applicable to all QCD regimes. Witten has 
generalized the approach to baryons [2!] and this has lead to a powerful l/N^ expansion 
method to study static properties of baryons, as for example, the masses, the magnetic 
moments, the axial currents, etc. The method is systematic and predictive. It is based 
on the discovery that, in the limit Nc 00, QCD possesses an exact contracted SU(2A^/) 



symmetry 



where Nf is the number of flavors. This symmetry is only approximate 



for finite Nc so that corrections have to be added in powers of l/Nc- The l/iV^expansion 



'jv^expan 



method has extensively and successfully been applied to ground state baryons [5|, [g, l7|, [Sj (for 
recent developments see Ref. js]). Its apphcability to excited states is a subject of current 
investigations. In this case the symmetry under consideration is assumed to be S\J{2Nf) 
X 0(3) where SU(2A^/) is related QCD, as introduced above. However 0(3) is not related 
to QCD but it brings an additional degree of freedom. It is of common practice to intro- 
duce it in order to construct orbitally excited states. The direct product SU(2A^j) x 0(3) 
is also used in quark models to classify three quark states, but there S\J{2Nf) is not an 
intrinsic symmetry. Thus the two approaches have formally the same symmetry in common 
which does not imply common dynamical assumptions. The only common feature is that 
the excited states are stable in a first approximation. 

The purpose of the present study is to see whether or not there is a compatibility between 
the two methods. If such a compatibility exists, an important support to the constituent 
quark model can be provided by the model independent 1/Nc expansion method, and a 
better understanding of the physical content of large Nc mass formulas can be gained. 
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In the language of quark models, the baryon states can roughly be classified into excitation 
bands with N = ior the ground state band and N = 1,2,3, .. . for excited states, where 
N represents units of excitation, like in a harmonic oscillator picture. The key tool of this 
comparative study is that one can analyze both the 1/Nc expansion results and the quark 
model basic ingredients in terms of N which makes the comparison between the two methods 
possible and very convenient. 

The paper is organized as follows. The next section introduces the mass formula used in 
the 1/Nc expansion method. Section III gives a mass formula obtained from a Hamiltonian 
quark model where the kinetic energy is relativistic, the confinement is an Y-j unction fiux 
tubes and the hyperfine interaction is of an one-gluon exchange nature. Section IV is de- 
voted to the comparison between terms of the mass formula which are common in the two 
approaches. The last section is devoted to conclusions. 

II. BARYONS IN LARGE A^^ QCD 

For simplicity, we illustrate the method with the Nj = 2 case but the arguments are 
similar to any Nf. So, here we deal with SU(4) which has 15 generators, the spin subgroup 
generators Si {i — 1, 2, 3), the isospin subgroup generators {a — 1, 2, 3) and Gia which act 
both on spin and isospin degrees of freedom. The SU(4) generators are components of an 
irreducible tensor operator which transforms according to the adjoint representation [211] 
of dimension 15 of SU(4). The SU(4) algebra is 

['S'i, Tg] — 0, [Si, Gja] — iSijkCrkaj [Ta, G^,] — i£abcGric, 

[Si, Sj] = iSijkSk, [Ta^Tfj] = iSabcTci 

[Giai Gjb\ = jSijEabcTc + J^ab^ijkSk- (l) 

Together with the generators ii of SO (3), the SU(4) generators form the building blocks of 
the mass operator. Then in the 1/A^c expansion the mass operator M has the general form 

M = J2ciOi, (2) 

i 

where the coefficients q are reduced matrix elements that encode the QCD dynamics and 
are determined from a fit to the existing data, and the operators Oi are 0(3) scalars of the 
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form 

c 

where Of^ is a fc-rank tensor in 0(3) and O'^^p a fc-rank tensor in SU(2)-spin (homomorphic 
to S0(3)), but invariant in SU(2)-flavor. Generally the operators are combinations of 
the SU(2A^/) generators and here, in particular, of SU(4) generators. The lower index i in 
the left hand side represents a specific combination. Each n-body operator is multiplied by 
an explicit factor of 1/A^"~^ resulting from the power counting rules j^. For the ground 
state, one has /c = 0. For excited states the k = 2 tensor is also important. The sum in 
the mass operator is finite. Operator reduction rules simplify the expansion. In addition, 
in practical applications, it is customary to include terms up to 1/Nc and drop higher order 
corrections of order As an example, in Eqs. (HI), we exhibit the list of operators used 

in the calculation of the masses of the [70,1^] multiplet up to order l/Nc included |lO |. 
Note that although O5 and Og carry a factor of their matrix elements are of order 

1/Nc because they contain the coherent operator G*" which brings an extra factor of A^^^. 

Oi = N,l, 02 = -^tS', O3 = ^T^T^ O4 = ^S'S\ 
' 

O5 = ii£(2)iiGmGia^ Oe = T'^G". (4) 

■'■^ c c 

Here Oi = N^t is the trivial operator, proportional to Nc and the only one which survives 
when Nc 00 [2], where the SU(4) symmetry is exact. It is the only spin-isospin inde- 
pendent term in the mass formula. The SU(4) quadratic operators S^S^, T'^T"' and G^^G*" 
should all enter the mass formula (the sum over repeated indices is implicit). But they are 
related to each other by the operator identity j3] 

{S\ S'] + {T", T"} + 4 {G'\ G'"] = ^Nc{3Nc + 4), (5) 

so one can express G^^'G*" in terms of S^S'^ and T°-T'^. Note that the right hand side of 
Eq. ^ is the eigenvalue of the Casimir operator for the irreducible representation [N^. — 1,1] 
of SU(4). The operators O2, O5 and Og are relevant for orbitally excited states. Among 
them, the role of O2 will be discussed below. 
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A. The ground state band 



The mass formula for the ground state up to order 1/Ncis simple because one can replace 
rparpa s^S^, due to an identity which holds for symmetric [Nc] states |7|]. As there is no 
orbital excitation, the mass formula ([2]) takes the following simple form 

M = c,iV. + c4^5^ + o(^), (6) 

which means that for = only the operators Oi and O4 contribute to the mass. Thus 
the fit gives quantitative information only for Ci and C4. For A^^ = 3, M^- = 940 MeV for 
S = 1/2, and Ma = 1232 MeV for S = 3/2, one gets 

ci = 289 MeV, C4 = 292 MeV. (7) 



B. Excited states 



Among the excited states, those belonging to the = 1 band, or equivalently to the 



[70Jj]_mu 

id 



13, 



18 



tiplet, have been most extensively studied, either for Nj = 2 jll 



12 



13 



14 



19 



20l | or for Nf = 3 [21]. In the latter case, first order corrections in SU(3) 



symmetry breaking were also included. 

The = 2 band contains the [56', 0+], [56, 2+], [70, i+] {i = 0, 2) and [20, 1+] multiplets. 
There are no physical resonances associated to [20, 1^]. The few studies related to the N = 2 



3and concern the [56', 0+] for iV. = 2 ^2^] , [56, 2+] for A^j = 3 |23j, andl70, i+] for Nf = 2 



2J] , later extended to Nf = 3 [25| . The method has also been applied [26| to highly excited 
non-strange and strange baryons belonging to [56,4"'"], the lowest multiplet of the A^ = 4 



band [27 1. 

The group theoretical similarity of excited symmetric states to the ground state makes 



the analysis of these states simple 



23 



261]. For mixed symmetric states, the situation is 



more complex. There is a standard procedure which reduces the study of mixed symmetric 
states to that of symmetric states. This is achieved by the decoupling of the baryon into an 



excited quark and a s ym metric core of A^r — 1 
the [70, l-lmultiplet (u 



12 



15 



16 



17 



quarks. This procedure has been applied to 



18 



19 



3, 



2l| and to the [70, i+] [i = 0, 2) 



multiplet 2J, |25| . In fact the decoupling is not necessary, provided one knows the matrix 
elements of the S\J{2Nf) generators between mixed symmetric states. The case of SU(4) 
has been presented in Ref. jlo| . 



In Section IV, we collect the values of Ci, C2 and C4 obtained in the above studies in order 
to make a comparison between those values and their analogs resulting from the quark model 
described below. 



III. QUARK MODEL FOR BARYONS 



A. Confining interaction 



In the framework of potential models, it is generally assumed that a baryon, viewed as 
a bound state of three quarks, can be described in a first approximation by the following 
spinless Salpeter Hamiltonian 

3 ; 

H = Y.^Jpf + m} + Vy, (8) 

i=l 

where mj is the current mass of the quark and Vy the confining interaction potential. 

The nonperturbative part of the gluon exchanges, responsible for the confinement, can 
be successfully described in the flux tube model 281]. In this framework, each quark is 
assumed to generate a string, or a flux tube, characterized by its energy density (string 
tension). Recent developments in lattice QCD tend to confirm the Y-junction as the correct 
configuration for the flux tubes in baryons [2^. In this picture, a flux tube starts from each 
quark and the tubes meet at the Toricelli point of the triangle formed by the three quarks. 
This point, denoted by xt, is such that it minimizes the sum of the flux tube lengths, and 
its position is a complicated function of the quark coordinates Xj. Moreover, the energy 
density of the tubes appears to be equal for mesons and baryons. The Y-junction potential 
reads 

3 

Vy = a \xi — xt\ ■ (9) 
1=1 



In Ref. [30[, it has been shown that this complicated potential is successfully approximated 



by the more easily computable expression 



V 



Oi ^^^^^ I iC 1 

i=l 



R 



'1 — a 



i<j 



(10) 



where R is the position of the center of mass. If a = 1, Eq. fllOp is a simplified Y-junction, 
where the Toricelli point is replaced by the center of mass. If a = 0, this interaction 



reduces to a A-type potential. Results of Ref. [30|, obtained in the framework of a potential 
model, show that a = 1 gives a better description than a = 0, and that the Y-junction is 
approximated at best by a close to 1/2. 



B. Mass formula 



Let us now introduce auxiliary fields, in order to get rid of the square roots appearing in 
the Hamiltonian (IHl). We get 



E 



2u, 



2 ^ 



2ui 



+ 



:i-a) 



E 

j<k 



a (X 



3 



2A 



+ 



A. 



jk 



•jk 



The auxiliary fields, denoted and Xij are, strictly speaking, operators. Although 

being formally simpler, H{fii, uj, Xij) is equivalent to H up to the elimination of the auxiliary 
fields thanks to the constraints 



S^^H{fii,Uj,Xij) = 
S^^H{fXi,Uj,Xij) = 
SxijH{ni,Uj,Xij) = 



^i,0 = 
\j,0 



a\xj — R\ 



(12a) 
(12b) 
(12c) 



It is worth mentioning that (/ij,o) can be seen as a dynamical mass of a quark of current 
mass mj, while (^'j,o) is^ in this case, the static energy of the straight string linking the 
quark i to the Toricelli point 3l||- Similarly, (Xij^) can be interpreted as the static energy 
of a straight string joining the quarks i and j. Although the auxiliary fields are operators, 
the calculations are considerably simplified if one considers them as real numbers. They 



are then finally eliminated by a minimization of the masses with respect to them [32]. The 
extremal values of /i,, uj, and Xij, considered as numbers, are logically close to the values 
of (/ij,o)) {^j,o)^ i^ijfi) given by relations ffT^ . This procedure leads to a spectrum 
which is an upper bound of the "true spectrum" (computed without auxiliary fields) 33| : 
it can be shown that, the more auxiliary fields are introduced, the higher are the masses 
compared to those without auxiliary fields 3^. Let us finally mention that, for a = 1, the 
Hamiltonian fllip can be related to the rotating string model for a baryon (see for example 



Ref. 



35|) 



In Ref . 36|] , it has been shown that the eigenvalues of a Hamihonian of the form fllip can 
be analytically found by making an appropriate change of variables, the quark coordinates 
Xi = {xi,X2,X3} being replaced by new coordinates = r/|. The center of mass is 

defined as 



R 



fllXi + ^2X2 + /i3^3 



(13) 



with /ii = /ii + /i2 + and ff} being the two relative coordinates. From Ref. [36|, it can 
be immediately found that the mass spectrum of bound states of three massless particles 
{rrii = for the u and d quarks) is given by 



3 / d _ 

M(/i, u, A) = uj{2n + £ + 3) + -ifi + au+ 



a] 



-A 



with 



Lu = a\l — 



a 3(1 — a) 
- + 

V 



2A 



(14) 



(15) 



n = + and i = £g + An obvious symmetry argument helps us to make the 
identification yUj = Vi = u, and Ajj = A. In this symmetric case, properties of the 
equilateral triangle together with the relations (fT2l) allow to make the following ansatz 



A = VSu. 



(16) 



Defining 



and 



we find 



Q = a 



Vs, v = Qu, 



N = 2n + 



M{fi, v) = aQ.fX{N + 3) + ^ (/i + z>) . 
\ Hiy 2 

Formula ( fT9l) is clearly symmetric in fi and z>. That means that we can set fi 
equality can be viewed as a sort of virial theorem. Then we have 

M(/i) = ^(iV + 3) + 3/i. 

One can easily find that the relation 5^M(/i) = implies 



(17) 

(18) 

(19) 
z>. This 

(20) 



1^0 



1/2 



(21) 



8 



and M(/io) = 6/io, as observed in Ref. [37|. Writing explicitly the square mass, we see that 
the model of Ref. [s?! also predicts Regge trajectories, which are in agreement with the 
experimental data for light baryons 

M^ifxo) = l2aQ{N + 3). (22) 

The Regge slope is here given by 12aQ. However, from experiment we know that the Regge 
slope for light baryons and light mesons are approximately equal. For light mesons, the 
exact value in the relativistic flux tube model is 27ra, a lower value than the one obtained 
from formula (!22l) . This is due to the auxiliary fields method: the more auxiliary fields 



we introduce, the more the masses are overestimated [S^]- What can be done to cure this 
problem is to rescale a: let us define a such that 12aQ = 27ia. Then, formula (l22ll is able 
to reproduce the light baryon Regge slope for a physical value a ^ 0.2 GeV^. Note that the 
best value for a is 1/2. Consequently, the best value for Q is 1/2 + v^/4 ^ 0.93. It is worth 
mentioning that such a rescaling of the string tension has already given good results in the 
study of hybrid mesons ^] . 

C. One gluon exchange and quark self-energy 

Although including only the confining energy is sufficient to understand the Regge tra- 
jectories of light baryons, it is well-known that the absolute value of the masses which are 
obtained are too high with respect to the experimental data. Other contributions are needed 
to decrease these masses and we shall estimate their effect perturbatively. The most widely 
used is a Coulomb interaction term of the form 

AM„,e = -^a.5^/-^\, (23) 

i<j * 

arising from one gluon exchange processes, where is the strong coupling constant, usually 



assumed to be around 0.4 for light hadrons [39[. Lattice QCD calculations also support this 



value [40l]. Assuming that {1/A) ^ 1/ {A), and using symmetry arguments, relations ( fT2|l 
lead to 

y^/ 1 \_3a^ V3aQ 
^ \ Ifj - Xj\/ Ao Ho 



AMoge = -2as^S-. (25) 
v3/io 



Another interesting contribution to the mass, which can be added perturbatively, is the 
quark-self energy. Recently, it was shown that the quark self-energy, which is created by the 
color magnetic moment of a quark propagating through the vacuum background field, adds 
a negative constant to the hadron masses |41ii| . Its negative sign is due to the paramagnetic 
nature of the particular mechanism at work in this case. The quark self-energy contribution 
for three massless quarks is given by 41 1 

AM,,e = -TT^. (26) 
27r/io 

The factor / has been computed in lattice QCD studies. First quenched calculations gave 



/ 



42]. A more recent unquenched work 



431] gives / = 3. Since its value is still a matter 



of research, we will only assume that / G [3,4]. 

With the unperturbed baryon mass M(/io), given by Eq. ( l22l) . the total mass is given by 
the sum Mq = M(/io) + AMoge + AM^^e- Then, in the first order of perturbation and for 
a = 1/2, it is straightforward to obtain the following mass formula for baryons 

4 12 
= 27ra{N + 3) - -^naa^ - ^^^^^ /a, (27) 

where the scaling 12aQ = 27icr has been used. The effects of the one gluon exchange term 
and of the quark self-energy are thus to shift the square mass spectrum by a global negative 
amount. Let us note that the symbol defined by Eq. (1181) and the quantity A^ used to 
classify baryon states and used to plot results from the 1/Ac expansion are the same. This 
common A^ will be used in the next section to perform a comparison between the results 
obtained in both approaches. 

The mass formula (1271) does not take into account spin relativistic contributions, as the 
spin-spin or spin-orbit forces. Within the auxiliary field formalism, all these corrections to 
';he static potential are expanded in powers of l//i^ where /i is the constituent quark mass 



441 ]. All the spin corrections to the mass formula fl27p must depend both on the matrix 
elements of the interaction and on the coefficient l//i^. In the following, we shall consider 
that the dominant dynamical effect is due to the constituent mass, while the matrix elements 
remain roughly constant with A^, as presented in the next section. 
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IV. COMPARISON OF THE TWO APPROACHES 



In the 1/Nc expansion method, the first term CiNc in the mass formula of Eq. ^ con- 
tains the main spin-independent contribution to the baryon mass, which in a quark model 
language, represents the confinement and the kinetic energy. So, it is natural to identify this 
term with the mass given by the formula (1271) . Then, for N^. = 3, we assume the relation 



e? ^ f . (28) 



which gives 



2% 

cl = —aN + co (29) 
y 

Fig. [1] shows a comparison between the values of c\ obtained in the l/N^. expansion method 
and those derived from the Eq. fl29p for various values of A^. From this comparison one can 
see that the results of large N^. QCD are entirely compatible with the formula fl2I?|) . From 
a fit, one has a = 0.163 ± 0.004 GeV^, a rather low but still acceptable value according to 
usual potential models, and cq = 0.085 ± 0.007 GeV^. To reproduce cq, we can set a<j = 0.4, 
/ = 3.5: these are very standard values. 

In most of the quark models however, the string tension is generally assumed to lie in the 
range [0.17, 0.20] GeV^. If the value of a is chosen in this interval, the corresponding values 
for cf, given by Eq. (1291) . are located in the shaded area of Fig. [T] Although the agreement 
with large Nc data is not so good than in the optimal case, where a = 0.163 GeV^, it remains 
satisfactory if we choose / = 3.98 (4.42) for a = 0.17 (0.20) GeV^, together with a, = 0.4. 
These values are larger than what is expected. It could be argued that other mechanisms 
than the quark self-energy are present, their contribution decreasing the total mass Mq. In 
mesons for example, retardation effects due to the finite interaction speed were shown to be 
also proportional to yU~^, like the quark self-energy 45|]. It is possible that, when retardation 
effects are included, / can again be chosen in the interval [3, 4] with a standard value of a. 
But, no model for retardation effects in baryon has been proposed yet. 

Within the auxiliary field formalism, we can expect that C2 and C4 oc /Xq ^, and thus 

^0 „o 
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FIG. 1: Values of cf computec 



(Eq. ^ ioi N = 0, Refs. 



21 



in the 1/Nc expansion (full circles) from a fit to experimental data 



261] for = 4), 



231] for iV = 1, Ref. [24] for iV = 2 and Ref. 
compared with results from a fit (see text) of the formula ()29p (empty circles and dotted line to 
guide the eyes). No data is available for iV = 3 in large Nc studies. Values of cf as predicted by 
formula (f^^ for a £ [0.17,0.20] GeV^ are located in the shaded area. 



We see that this behavior is coherent with the large Nc results in Figs. [2] and [3l We chose 
= 208 ± 60 MeV so that the point with = 1, for which the uncertainty is minimal, is 
exactly reproduced. Let us note that the spin-orbit term is vanishing for = 0, so no large 
Nc result is available in this case. To compute the parameter a fit is performed on all the 
large A^^ data. We obtain then = 1062 ± 198 MeV. Note that > c^. This shows that 
the spin-spin contribution is much larger than the spin-orbit contribution, which justifies 
the neglect of the spin-orbit one in quark model studies. 



V. CONCLUSIONS 



This study supports the quark model basic assumptions by the compatibility of its mass 
formula with the mass formula derived from the model independent 1/Nc expansion. These 
assumptions are: relativistic kinetic energy for light quarks, Y-junction confining interaction, 
negligible spin-orbit interaction, hyperfine interaction dominated by a spin-spin term. A 
recent analysis shows that a fiux tube model and a feeble spin-orbit interaction give a 
succesful account of hadron spectroscopy 
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FIG. 2: Values of C2 computed in the \/Nc expansion (full circles) from a fit to experimental data 



(Refs. [2l|,l2^ for = 1, Ref. 



2]] 



for N = 2 and Ref. [26] for = 4), compared with results from 



formula (j3ip (empty circles and dotted line to guide the eyes). No data is available for = 3 in 
large N^. studies. 

In addition this study suggests that a good description of the bulk content of the baryon 
mass can be obtained with a spin independent energy eigenvalue of the form Mq oc a/A?" + 3 
where A^ = 0, 1,2, . . . is the number of excitation units, as in the harmonic oscillator. It 
also shows that the spin-orbit and spin-spin interactions vanish with the excitation energy. 
Moreover this comparative study gives a better insight into the large Nc mass operator where 
the coefficients Cj encode the QCD dynamics. 
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